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,; Abstract. The variation of an ionizing flux as a mechanism to stimulate the condensation of a diffuse gas is considered. 
' To illustrate this effect, two situations are examined: one on the context of pregalactic conditions, and the other on 
, the context of the actual interstellar medium. We focus our attention on flash-like variations; that is, during a "short" 
period of time the ionizing flux is enhanced in comparison to the pre- and post-flash values. In both cases the cause 
of the induced phase change is the same: the enhancement of the cooling rate by the increase in the electron density 
caused by the momentary increase of ionizing flux. After the passing of the flash, the cooling rate remains enhanced 
^f) , due to the "inertia of the ionization" . In the first case (metal free gas) the cooling rate is enhanced due to the fact that 
the increase of the electron density makes possible the gas phase formation of H2 by the creation of the intermediaries 
and i?^ . We show that after the passing of the photo-ionizing flash a cloud near thermo-chemical equilibrium at 
~ 8000 K may be induced to increase its H2 content by many orders of magnitude, causing a rapid decrease of its 
temperature to values as low as 100 K. In the second case (solar abundances gas) the dominant cooling mechanism of 
' the warm neutral gas (the excitation of heavy ions by electron impacts) is proportional to the electron density. We 
show that, for the expected states of the warm interstellar gas, ionizing flashes may induce the phase transition from 
O the warm to the cool phase. The results indicate that the mechanism of induced condensation studied here might 
^ I play a relevant role in the gas evolution of the diffuse gas in both, the pregalactic and the actual interstellar medium 
^ ■ conditions. 
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1. Introduction 

Thermal condensation of diffuse gas is a commonly invoked route for the formation of dense astrophysical structures. 
This thermal condensation may occur spontaneously or may be induced by variations of the external conditions. In the 
first situation (spontaneous condensation), the initially diffuse gas is in a state out of equilibrium in which the cooling 
dominates over the heating. Eventually, the gas reachs a thermo-chemical equilibrium in a dense and cool state. In 
the second situation (induced condensation) , in absence of variation of external conditions the gas would reach (or be 
in) a state of thermo-chemical equilibrium in a diffuse and hot state, but in presence of appropriate variations of the 
external conditions, the gas evolves toward a cooler and denser state. Induced condensation of diffuse gas must play 
an important role in the large scale evolution of the medium; specially if the stimulating sources are able to induced 
condensation far away from them. The synchronization and the large scale patterns of star formation may be governed 
by this kind of stimulation. In fact, most of the models and numerical simulations of the star formation in disk galaxies 
include as a main process the so-called self-propagating star formation (Gerola & Seiden 1978; Seiden & Gerola 1982; 
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Shore 1981, 1983; Palous et al. 1990; Comins & Shore 1990; Cammerer & Shchekinov 1994). This kind of systems are 

part of the wider class of reaction-diffusion systems (Kapral 1993, and references therein). In particular, the so-called 
exitable media are appropriate to capture the main features of the pattern formation in disk galaxies (Smolin 1996). 

Among the external conditions whose variations would induce phase transitions are the external pressure and the 
ionizing flux. We focus our attention on the effect of the variation of the ionizing flux, keeping in mind that induced 
condensation by pressure variations may play a relevant role (Roberts 1969; Shapiro & Kang 1987). In particular, we 
focus our attention on flash-like variations; that is, during a "short" period of time, the ionizing flux is enhanced in 
comparison to the pre- and post-flash values. To illustrate this mechanism of induced condensation, two examples are 
given here: one on the context of pregalactic conditions (in a free metal cloud), and the other on the context of the 
actual interstellar medium (in a gas with solar abundances). In both cases the cause of the induced phase change is 
the same: the increase in the electron density due to the momentary increase of ionizing flux enhances the cooling 
rate. After the passing of the flash, the cooling rate remains enhanced due to the inertia of the ionization (i.e. the 
characteristic recombination time is much larger than the cooling one). 

Many scenarios in which a flash of radiation affects the evolution of gas clouds can be imagined. For the metal free gas 
case, this kind of induced condensation might be relevant at the epoch of galaxy formation. In particular, the ionizing 
flash effect might be a key step in the sequence of events that have conduced to the formation of globular clusters. It has 
been stated by Cox (1985) that "such a rapid fragmentation of the halo almost certainly requires inducement by energy 
leaving the (primitive) disk" . For the case of a gas with solar abundances, the evoked mechanisms of induction are 
diverse but generally associated to compression and convergent mass flows (see review by Elmegreen 1992). However, 
in addition to these mechanisms, the induced condensation by ionizing flashes appears as an effect to be considered. 

In Sec. 2 the basic equations to follow the thermal and chemical evolution of a gas subject to variations of the gas 
pressure and the ionizing flux are given. In Sec. 3 the evolution of a free metal gas cloud subject to a flash of ionizing 
and dissociating radiation is analyzed, whereas, in Sec. 4, the evolution of a solar abundance gas subject to a flash of 
cosmic ray flux is considered. Finally, the conclusions are given in Sec. 5. 



2. Basic equations 

The thermal and chemical evolution of a gas subjected to variable external conditions are calculated by solving the 
system of kinetic equations and the equation of energy conservation. The system of kinetic equations can be written 
in the general form 

(1) 

+ (^mXm - Q-iXi , 

where n is the total number density of hydrogen nuclei and Xj is the relative number density of the z-th species. In 
Eq. (1), ajk and an are, respectively, the formation and destruction rates of the z-th species due to double collisions, 
whereas, and are, respectively, the production and destruction rates of the i-th species due to the interaction 
of particles with radiation. 

The equation of energy conservation equation can be expressed as 

where U is the internal energy by unit mass, C is the net cooling rate by unit mass (C ~ {K — V)/p), P is the pressure, 
p is the mass density, and V is the speciflc volume. Assuming that P = pRT/ ji and U = 3RT/2iJ,, with R the gas 
constant, and /j. the molar mass, eq. (2) can be written as 

dT^Tdp 2TdP_2ii 

dt iJ. dt 5 P dt 5 R ' ^ ' 

Two extreme situations that can be used to conflne intermediate situations are the constant density and constant 
pressure approximations. In these approximations, the energy conservation equation can be expressed as 

^ = ^^-A^£, (4) 
dt fj, dt R ' ^ ' 



where, A = 2/3or2/5in the constant density or in the constant pressure approximations, respectively. 
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Variations of the external ionizing sources flux produces direct variation on the rates am and ai, and on the heating 

rate r, whereas, variations of the prcsurc affect density and temperature. Notice that variations of n or ai provoke a 
variation of the chemical state of the gas, and therefore produce a variation of the cooling rate A. In the following, 
we solve the basic equations (1) and (2) with the appropriate assumptions in order to model: a) the evolution of a 
metal free cloud subjected to a variation of the ionizing and dissociating flux (sec. 3), and b) the evolution of a solar 
abundance gas subjected to variations of the primary ionization rate due to cosmic rays (sec. 4). 

3. Photo-ionizing radiation flash as trigger of efficient cooling in free metal gas clouds 

In the context of formation of galaxies and globular clusters, the non-equilibrium formation of H2 has been identified 

as a key process for achieving a rapid cooling bellow lO^K (Palla & Stabler 1983; Izotov & Kolcsnik 1984; Shapiro & 
Kang 1987; Palla & Zinnecker 1987; Kang et al. 1990; Anninos et al. 1996; Padoan et al. 1996; Tegmark et al. 1996). 
On the other hand, the ionization and dissociation of primordial gas by UV background radiation have been taken into 
account as an important parameter in many studies of the thcrmo-chcmical evolution of pregalactic and intcrgalactic 
structures (Kang et al. 1990; Donahue & ShuU 1991; Ferrara & Giallongo 1996; Haardt & Madau 1996; Navarro & 
Steinmetz 1996; Mucket & Kates 1997). In general, the UV background radiation acts as an inhibitor of H2 formation. 
However, there are exceptions to this rule. Recently Haiman ct al. (1996) have shown that UV background radiation 
can enhance the formation oi H2 in primordial gas at high densities (> Icm""^) and low temperatures (< lO^K); but for 
densities lower than the above value, the effect of a constant UV background radiation is to inhibit the H2 formation. 
As it will be shown in the present section, rapid variations of the UV background radiation can also enhance the 
non-equilibrium formation of i?2- In particular, we focus our attention on flash like variations capable of heating and 
increasing the ionization fraction of an initially warm-neutral gas cloud near thermo-chemical equilibrium. After the 
passing of the radiation pulse, the enhanced ion fraction makes the gas phase formation of H2 molecules possible by 
the creation of the intermediaries H~ and H^- The presence of small quantities of H2 molecules then makes possible 
further radiative cooling to temperatures as low as ~ 10^ K. The rapid cooling (^ 0.2 free-fall times) abruptly reduces 
the Jeans mass by a factor ^ 10^, permitting the fragmentation of clouds initially marginally stable. Izotov (1989) 
has considered the homogeneous contraction approach of a gravitational unstable cloud, but here we are interested 
in delimiting the necessary conditions that provoke the rapid cooling of the cloud without invoking the gravitational 
collapse. UV radiation pulses have similar effects than shock waves because the post-shock flow also recombines out 
of equilibrium. 

The spectrum, amplitude and duration of the radiation pulse are free parameters in our model. However, it should 

be noticed that, in order to induce efficient H2 cooling in the cloud, the detailed form of the flash is not important 
as long as that, during the flash, the ionization of the cloud increases appreciably and after this, the ionizing flux 
decreases to background values in a short time compared to the recombination time. 

The initial conditions of the cloud and the characteristics of the hypothetic ionizing pulse depend on the chosen 
scenery. For example in the Fall and Rees (1985) scenery, during the proto-galactic collapse, if the gas is assumed to be 
lumpy, the overdense regions will cool more rapidly than the underdense regions producing a two-phase medium. But 
the developing rate of this overdense regions depends on the initial conditions (i.e. the denser regions develop faster; 
Murray and Lin 1990) and on their interaction with the system. Thus, a dependence of the density contrast and of the 
mean cloud masses on the galactocentric distance are expected. The assumption of an ionizing flash produced in the 
galactic nucleus implicitly assumes a radial increment of the delay in the evolution of the clouds relative to the center. 
The dilution and attenuation of the ionizing and dissociating radiation produced in the flash also introduce a radial 
dependence. In fact, the possibility that proto-galactic structures are exposed to UV radiation emitted by massive 
young stars or an active galactic nucleus has been considered previously (Kang et al. 1990). Therefore, it would be 
useful to know the dependence of the cloud evolution on the flash characteristics, and on the cloud initial state (i.e 
cloud mass (Mc), temperature (T), number density (n), and the relative number density of the species (Xj). 

To follow the thermal and chemical evolution of a metal-free cloud in presence of a variable radiation flux, we 
consider an idealized uniform cloud. The gas model adopted for this application considers the following 9 species: H, 
e, H~, H2, H2, He, H^, and H^~^. Assuming a mass fraction y (= 0.3) of Helium relative to Hydrogen, then 
H = {l + y)/{l + Xe — XH2 ^-^H2 +y/3-9'''), where Xi is the ratio of the number density of species i to the total number 
density of Hydrogen nuclei n. The processes of formation and destruction of these nine species are assumed to be the 24 
reactions in Table 1 of Rosenzweig et al. 1994. The adopted cooling rates, include: 1) the collisional ionization of H by 
electron impact; 2) the free-free transitions of (Izotov (1989) and references therein); 3) the collisional ionization 
by electron impact of and if+; 4) the free- free transitions of and H^^; 5) the total dielectronic cooling rate of 
(Shapiro & Kang 1987); 6) the H2 cooling rates due to rotational and vibrational transitions excited by [H — H2) 
and {H2 — H2) collisions, calculated according to Lepp and Shull (1983). These H2 cooling rates must be multiplied 
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by a factor L{Tvib) and L{Trot) to take into account the escape probability of the vibrational and rotational photons, 
respectively. Finally, 7) the excitation of the H2 vibrational levels {v < 4) by low energy electron-impact was taken 
into account (Klonover & Kaldor, 1979). 

The external radiation flux affects the photoionization of the Hydrogen and Helium, and the H2 photo-dissociation 
(and its associated heating rates; see Rosenzweig et al. 1994). The external photon flux is assumed to have the 
quasar-like distribution (Shapiro and Kang 1987) 



where uh is the H Lyman-edge frequency. The function e{t) is introduced to represent the changes of the flux level 
outside the cloud, and is assumed to behave as eion(t) in the Lyman continuum, and as edis{f^ for the photo-dissociating 
photons. Since a homogeneous cloud is assumed in this application, the optical depth at its center, for the ionizing 
radiation, can be written as: 



where M is the mass of the cloud in units of and al are the absorption cross sections of species i at frequency 
V. As a first approximation, the attenuation of the external radiation field t) by the factor exp""^" is adopted to 
schematize the opacity effects. For the dissociating radiation field a similar approximation is adopted (Rosenzweig et 
al. 1994). Obviously, by using this first approximation, the effect of the ionizing radiation is minimized, and the study 
of the internal structure of the cloud is not possible. 

For the time dependence of the flux level e we assume a flash like variation. More precisely, during a first period of 
integration the external radiation flux remains at the assumed background value (i.e. Cjon = ^dis = ^o)- The duration 
of this initial period of integration is large enough (i.e. ^ 10^ yr) to ensure that the cloud is near its thermo-chemical 
equilibrium state. After this first period of stabilization, it is assumed that both, Cjon and e^is, increase linearly to the 
value ei in a time Tincr, and remains at this value during a time t flash- Then, the flux level is assumed to decrease 
exponentially with a characteristic time Tdecr to the background value eo for the ionizing flux, and to £2 for the 
dissociating flux. This kind of flux variations are expected in the case when the radiation pulse is assumed to be 
produced by an intense but short event of stellar formation or by the radiation coming from a front shock. In a few 
million years after the end of the stellar formation process, the ionizing flux is expected to decrease to background 
values. On the other hand, due to the contribution of intermediate and low mass stars to the dissociating flux, it is 
expected that its level remain higher than the background value (e2 > eo) long after the death of massive stars. 

For a given initial state of the cloud, there is a critical background flux level ecri, bellow which the cloud spon- 
taneously cools because the self- shielding allows the formation of H2 at relatively high concentrations. This critical 
background flux level is mainly sensible to the initial ion fraction due to its strong influence on the H2 formation rate. 
In order to describe this initial state dependence, it is useful to look at the thermo-chemical equilibrium (TCE) curve 
(i.e. £ = 0, and ^ = 0). Figure (l.a) shows the TCE curves for two different values of the background flux level 
(lower curve eo = 1) and upper curve eo = 10) when the mass of the cloud is lO'' M©. Note the three phase structure 
characterized by the presence of three stable branches (solid lines) denoted in Fig. (l.a) as W-I (i.e. warm-Ionized; 
T Z 13000K), W-N (i.e. warm-Neutral; 8000K <, T <, 6000K), and C-N (i.e. cool-Neutral; T <, 300K). Note also 
the qualitative difference of the TCE curve for low (eo = 1) and high (eo = 10) background flux level. That is, for 
eo = 10 the left-hand maximum (denoted as Pmax-i) is over the right-hand maximum (denoted as Pmax-2), whereas, 
Pmax-i < Pmax-2 for eo = 1. This difference is important because if the cloud is initially in TCE in the W-I branch, 
and the pressure is progressively increased from P < Pmax-i to P > Pmax-i a transition to the W-N phase occurs 
in the case eo = 1, but in the case eo = 10 the C-N phase is reached. For eo = 1, Fig. (l.b) shows the typical isobaric 
evolution tracks for initially ionized states (arrows A, B, C and D) and for initially neutral states (arrows E, F and 
G). For cases A, B, and D, as expected the cloud evolves toward the stable branch W-I, W-N, and C-N, respectively. 
However, in case C, even when P < Pmax-2 the track reaches the C-N branch transversing the W-N branch due to 
the inertia of the ionization fraction. More precisely, the excess of the ionization fraction shifts the maximum pressure 
Pmax-2 of the Corresponding thermal equilibrium curve to a value bellow the pressure corresponding to case C. On 
the other hand, for the initially neutral states the evolution is substantially different. In case E the track stops in the 
W-N branch; note the difference with case C. In case F the cooling rate is so small that in practice the cloud reaches 
a quasi-stationary state similar to that in the W-N branch. The small cooling efSciency is due to the very small H2 
formation rate for the low ion fraction in the quasi-stationary W-N state. Finally, for high enough gas pressure as in 




(5) 



= 6.6 X 10i8n(af + a^^Xn, + erf =^X^+)(— )V3 



(6) 
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Fig. l.a. The thermo-chemical equilibrium curves for a lO'^ Mq cloud exposed to the two labeled values of the background flux 
level (see text). 




log nH(cm 



Fig. l.b. The thermo-chemical equilibrium curve and the typical isobaric evolution tracks for initially ionized states (arrows 
A, B, C and D) and for initially neutral states (arrows E, F and G). 



case G, the cloud reaches the C-N branch. The isobaric evolution tracks for a high background flux level (i.e. eo = 10) 
are similar to those in Fig (l.b) except that the W-N branch can not be reached from an ionized initial state as in 
case B in Fig. (l.b). It is interesting to note that when the cloud is in the quasi-stationary state schematized by track 
F, a large enough ionizing flash is able to increase H2 formation rate and to stimulate the rapid condensation to reach 
the C-N branch. 

The analysis of the various evolution tracks in Fig. (l.b) have been made for a fixed background flux level. Taken 

into account that an increase of the background flux level shift the TCE curve upward, a similar analysis can be made 
for a fixed pressure but varying the background fiux level. Figure (2. a) shows, in the constant pressure approximation, 
the time dependence of the gas temperature for a 10^ Mq cloud subject to the three labeled values of the background 
flux level eo. The initial condition is the same in the three cases and corresponds to an out of equilibrium warm-ionized 
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Fig. 2. a. The time dependence of the temperature for a lO'^ Mq cloud exposed to the three labeled values of the background 
flux level and an initial out of equilibrium warm- ionized state (see text). 
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Fig. 2.b. The time dependence of the temperature for a 10 Mq cloud exposed to the two labeled values of the background 
flux level and an initial out of equilibrium warm-neutral state (see text). 

state at a gas pressure log{P/k) = 2.046, and n = 2.5 x IQ-^cm-^, T = 2 x lO^K, Xh = 8.8 x 10-*, Xh^ = 2.3 x lO'^^, 
Xh+ = 0.999, Xe = 1.147, Xh- = 4.6 x lO'^S X„+ = 1.2 x lO'^S Xhs = 2.57 x 10-^ XHe+ = 3.11 x 10-^ and 
XHe++ = 7.24 X 10~^. The temperature evolution in the three plotted cases differs because for eq — 1-6 the system 
stabilizes in the W-I branch (i.e. as in track A of Fig. (l.b)), for eo = 1.1 the system stabilizes in the W-N branch (i.e. 
as in track B of Fig. (l.b)), and for eo = 1.0 the system stabilizes in the C-N branch (i.e. as in track C of Fig. (l.b)). On 
the other hand, Fig. (2.b) shows the temperature evolution for the same pressure as in Fig. (2. a), but for an initially 
neutral warm state, (i.e. n = l.Ox lO'^cm'^, T = 1.03x lO'^K, Xjy = 0.9992, X^, = 5.801 x 10"^, Xh+ = 7.978x10-", 

Xe = 2.631 x 10-3, Xh- = 3.914 x 10-^, X„+ = 4.813 x 10-^ Xhc = 7.374 x 10-^, XHe+ = 1-819 x 10-^ and 

2 

XHe++ = 7.169 X 10-^). In this case, the evolution for eo = 0.3 corresponds to the situation schematized by track 
F in Fig. (l.b); that is, the system enters in a quasi-stationary state characterized by a very slow decrease of the 
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temperature. However, for eo = 0.2 the system evolves rapidly to the C-N branch, a situation that corresponds to that 

schematized by track G in Fig. (l.b). The results in Figs. (2. a) and (2.b) illustrate how the initial cloud state affects 
the critical value of eo bellow which the cloud evolves toward the C-N branch (i.e. the cloud condenses). In order to 
show the effect of the constant density approximation. Fig. (2.b) also shows the temperature evolution for eo = 0.2 
when the constant pressure approximation is switched to the constant density approximation at the time when the 
rate of decrement of the cloud radius equals the sound speed. 




For initially ionized states at pressure P, the critical value Cch roughly corresponds to the value of eo that produces 

a maximum in the TCE curve at pressure P. For a high background flux level, the left-hand maximum (denoted 
as Pmax-i in Fig- (l-a)) is the higher maximum and must be equaled to P. For a low background flux level, the 
right-hand maximum must be equaled to P. On the other hand, for initially neutral states, the critical value Cch must 
be calculated by finding the value of eo bellow which the evolution does not fall in the quasi-stationary warm-neutral 
state. The dependence of the critical value of ecri on the initial cloud state is summarized in Fig. (3) for a 10'' Mq 
cloud. The results are showed as function of the initial cloud column density [N = n x Rdoud)- The left-hand side 
curve corresponds to the critical background flux level for initially ionized states. The upper segment (labeled I-branch) 
corresponds to background flux levels for which Pmax-i > Pmax-2, and therefore, the condition Pmax-ii^cri) — P is 
used. The lower segment (labeled N-branch) corresponds to the case when the condition Pmax-2{^cri) = P is used 
because Pmax-2 > Pfnax-i for these values of Ccri- The initial column density is calculated assuming that initially the 
cloud is at 2.0 x lO^'K and at a pressure P = Pmax{^cri)] the initial concentrations are assumed to be the used in Fig. 
(2. a). For values of {N,eo) above the left-hand €cri{N) curve, the initially ionized cloud reaches TCE in the stable 
warm phase, but below this curve the cool-neutral branch is reached. On the other hand, the right-hand side curve 
corresponds to the critical background flux level for initially neutral states. In this case, the initial column density 
is calculated assuming that initially the cloud is at 1.2 x 10*K and the initial concentrations are assumed to be the 
used in Fig. (2.b); the initial density is varied in order to cover the plotted range of ecri- For values of (A^, eo) above 
the right-hand ecri{N) curve, the initially neutral cloud attains a quasi-stationary warm-neutral state, but below this 
curve the cool-Neutral branch is reached. It is to be noticed that as the initial electron concentration increase, the 
right-hand ecri{N) curve approaches the left-hand curve. In any case, there exist a set of initial conditions for which 
the cloud reaches a quasi-stationary warm-neutral state. The point to be emphasized is that these quasi-stationary 
states are susceptible to be induced to condensate if the cloud is exposed to an intense enough ionizing flash. 

In order to illustrate the effect of the ionizing flash on a cloud that has fallen in the quasi-stationary warm- 
neutral state, Figs. (4a,b,c) show respectively the time dependence of the temperature, the relative number density 
of electrons (-^e-)) and of Hydrogen molecules {XH2), for a IO^'Mq cloud. The initial condition corresponds to an 
out of equilibrium warm-neutral state at a gas pressure log{P/k) = 2.074 (i.e. n = 1.1 x 10~^cm~^, T = 1.0 x lO^K, 
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Fig. 4a. The effect of the ionizing flash e{i) on the the time dependence of the temperature for a 10^ M0 cloud. The curve 
labeled (a) corresponds to the case e(t) = eo > €cri- The curves labeled (b) and (c) correspond to the evolution of the cloud 
when it is subjected at t = 2 x 10* yr to a flash (see text). 




Fig. 4b. The time dependence of the relative number density of electrons {X^- ) for the initial and external conditions corre- 
sponding to Fig. (4a). 



and the initial concentrations used in Fig. (2. a)). The curves labeled (a) in Figs. (4) are plotted for reference, and 
correspond to the case when the cloud is subjet to a constant background flux level of eo = 1.0. As expected for 
e{t) = eo > Scri, the cloud evolves toward a quasi-cquilibrium warm neutral state. The curves labeled (b) and (c) in 
Figs. (4) correspond to the evolution of the cloud when it is subjected at i = 2 x 10® yr to a flash. In case (b), the 
flash characteristics are: Tincr = 5 x 10^ yr, Tjiash = 2 x 10^ yr, r^ecr = 2 x 10^ yr, ei = 10^, and €2 = 1.5 x eo. In 
case (c) the flash characteristics are the same as in case (b) but with e2 = 2 x eo. Notice in Figs. (4) that during the 
increase of the external radiation flux from eo to ei, the electron density and the temperature increase, whereas the 
H2 and H~ densities decrease. During the time when e{t) = ei, the electron density and the temperature continue 
to increase because the variation of e{t) during the time r,„cr is rapid enough to leave the gas far from equilibrium. 
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Fig. 4c. The time dependence of the Hydrogen molecules {XH2), for the initial and external conditions corresponding to Fig. 
(4a). 

At the end of the lapse of decrement of e{t), if Tdecr is short enough, the gas has an excess of electrons and thermal 
energy compared with the equilibrium values corresponding to Cjon = eo and eaiss = £2- After the lapse of UV flux 
decrement, recombination continues, but at a lower rate than the cooling. The excess of electrons at these relatively 
low temperatures results in an enhancement of the H2 rate formation. Even when e2 > cq, the abundance of H2 may 
reach a large enough value to produce considerable self-shielding. If the cloud reaches a critical value (r ^ 1) for 
the optical depth at dissociating frequencies, the H2 abundance grows very fast, allowing the cooling of the cloud to 
temperatures of the order of lOOK. This is the situation for case (b) in Figs. (4), where the flash induce rapid cooling 
even when £2 = 1.5 x cq. On the other hand, for case (c) a post flash dissociating level with e2 = 2 x cq is enough to 
inhibit the formation of i?2, and then, the cloud remains warm. 

A detailed study of the dependence of the cloud evolution on the phase space of free parameters (cloud mass, 
initial conditions, and flash characteristics) is out of the scope of this simple application. However, such detailed study 
may reveal that the stimulating condensation process studied here can be effective in a restricted region of the free 
parameters space, and therefore, may act as a selective effect that contributes to the formation of dense structures at 
certain scales. 



4. Ionizing pulse as trigger of warm to cool phase transition in a gas with solar abundances 

Among the various processes that determine the change of state of the ISM material, phase transitions are expected to 
play an important role. In particular, warm to cool phase transition has been indicated as a channel to transform diffuse 
gas (< lO^^cm^"*) into denser states (> lO-'^cm"^) (Field, Goldsmith & Habing 1969; Lepp ct al. 1985; Parravano 
1987; Lioure & Chieze 1990; Dickey & Brinks 1993). Moreover, it has been proposed (Parravano 1988, 1989) that the 
large scale star formation rate must be self-regulated because, on one hand the warm gas condensation is inhibited 
by high enough UV radiation (coming mainly from massive stars), and on the other hand, a gas supply from the 
diffuse phases is required to feed the large scale star formation process. In this way, the large scale star formation 
rate is limitated, and the warm gas tends to remains close the critical state for warm to cool phase transition. The 
fact that large quantities of ISM warm gas remain close to this critical state allows the trigger of condensation by 
relatively small variations of the ambient conditions. More precisely, the mean "distance" of the warm gas state to the 
critical state for spontaneous condensation is determined by the amplitude spectrum of the variations of the ambient 
conditions. In any case, the study of triggered condensation of warm gases close to the critical state can be justified 
by the self-regulatory hypothesis. It is to be noticed that the scale of the inhibitory process is expected to be much 
larger than the scale of the triggering processes: a fact that apparently is a common characteristic of many dynamical 
systems where spiral structures arise (Smolin 1996). 
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Triggering mechanisms of star formation are usually related to compression and pushing of the ISM gas by high 

pressure events associated to stars formed previously (reviews on this topic can be found in Elmcgrcen 1992, Franco 
1992). Also, compression of the warm gas by the spiral density wave have been evocated as a main trigger of its 
condensation (Roberts 1969); however, stimulated condensation by ionizing flashes may also enter as an initiator 
of the chain of processes that finally results in the formation of stars. As it will be shown bellow, this stimulating 
mechanism is particularly efficient when the warm gas state is close to its marginal state for spontaneous phase 
transition. Moreover, variations of the ionizing flux are expected to precede pressure variations if both variations are 
associated to the same perturbing event. Large local variation of the ionizing rate arc expected to be present in the 
interstellar medium. Sudden appearance and disappearance of ionizing sources occur continuously in the galactic plane. 
Also, "rapid" variation of the opacity to the ionizing radiation in a line of seeing is expected. Finally, variations of the 
cosmic ray flux are expected to be present due to variations of the sources (cosmic ray acceleration in shock fronts 
with oblique B-fields; Blandford and Ostriker 1978), and to variations of the magnetic field topology (i.e. focalization 
or dis-focalization of the cosmic ray stream in a region). It was also proposed (Ko & Parker 1989; Nozakura 1993) 
that star formation controls dynamo activities and hence large scale magnetic fields of disk galaxies. 

The cooling of the diffuse warm interstellar gas is mainly due to electron collisional excitation of a) fine structure 
levels and metastable states of the positive ions C+, Fe+, and 5+, and b) Lya excitation. Therefore, an increase 
of the ionizing flux (and the consequent increase of the electron density) tends to reduce the gas temperature if the 
increase of the cooling efficiency overcomes the associated increase of the heating. In general, this is the case for low 
equilibrium ionization fractions when the kinetic energy of the electrons that result from the ionization process is 
small. Depending on the value of the pressure and on the increase of the ionizing flux, the new thermal equilibrium 
state might be located in the cool neutral branch. That is; sufficiently large variations of the ionizing rate are expected 
to provoke the warm to cool phase transition. Moreover, if the ionizing flux variation has a flash-like variation, then, 
after the passing of the flash, the cooling rate remains enhanced due to the inertia of ionization. 

In order to illustrate this effect, we consider variations of the primary ionization rate by cosmic rays C,o which is 
the accepted main source of ionization of the interstellar warm gas far away from massive stars. Other mechanisms of 
ionization have been proposed (i.e. OB stars (Reinolds & Cox 1992), the neutrino decay theory (Sciama 1990, 1993)), 
but for the present analysis the exact mechanism of ionization is not relevant. What is important here is the variation 
of the ionizing and heating rates due to a change of the considered ionizing flux. 




2.2 - only Warm Phase 
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Fig. 5. The maximal {Pmax) and the minimal {Pmin) pressures of the thermo-chemical equilibrium curve plotted as a function 
of (see text). 

The main changes of the curve of thermo-chemical equilibrium (TCE) when C,^ is changed are summarized in Fig. 
(5), where the maximal {Pmax) and the minimal {Pmin) pressures of the TCE curve are plotted as a function of Co- All 
the results in this section correspond to the standard solar neighborhood far UV energy density and gas composition. 
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The labels are used to remind that if P < Pmin (or P > Pmax) then there is only one possible state of TCE in the 

diffuse and warm phase (or in the dense and cool phase). If Pmin < P < Pmax then the gas may reach any of the two 
stable branches (the warm or the cool branch). If the gas is initially in TCE at the warm branch (for example with 
the external conditions (PqjCoa) corresponding to the point (A) in Fig. (5), then, the gas could be forced to evolve 
toward the cool branch if the external conditions are changed, for example, to (PotCob) corresponding to the point 
(B). The typical time for the transition from the warm to the cool branch is 10^ yr. Once the gas reaches the cool 
phase, the external conditions can change again to (Pq, Coa) but the gas will remain in the cool branch. To drive the 
gas to the warm phase again, the external conditions must be changed, for example, to (PoCiCoc) corresponding to 
the point (C) in Fig. (5). Now, the gas could return to the initial state if the external conditions changes to the initial 
condition (Po,Coa)- 

The phase transitions described above assume that the time between the consecutive changes of Co are long enough 
to reach TCE. If the variation of Co occurs before TCE is achieved, then, the phase transition does not necessarily 
occur. Here we will consider the effect of flash-like variations of Co at constant pressure. That is, at the beginning we 
assume that the gas is in the warm branch in the equilibrium state corresponding to the external conditions (Pq, Coa)- 
Then, Co is increased abruptly by a factor during a lapse of time Aj, after which the primary ionization rate by 
cosmic rays returns to the initial value Coo- As mentioned above, a fact that favors the phase transition is that the 
recombination time is much larger than the cooling time; then after the end of an ionizing pulse, the cooling rate 
remains enhanced due to the non-equilibrium excess of electrons. 

In this simple analysis we neglect non-local processes as thermal conduction, cosmic ray attenuation, radiative 
transfer, and gas dynamics. Only local processes of radiative cooling, heating, H ionization and recombination are 
considered (Parravano 1986), assuming that the pressure remains constant during evolution. The neglected non-local 
processes may play an important role in the local evolution of the gas and in should determine the spatial variations 
of the physical variables. However, here we are interestested in showing that the inertia of ionization (after the passing 
of the ionizing pulse) can, in many cases, enhance sufHciently the cooling efficiency to produce a phase transition from 
the warm gas phase to the cool phase. The evolution of the temperature and the ionization degree is calculated by 
solving simultaneously the energy conservation equation (2), and the kinetic equation (1) restricted to the Hydrogen 
ionization-recombination processes. 




Fig. 6. The temperature evolution of the warm gas initially in thermo-chemical equilibrium close to the critical state for 
spontaneous condensation. The three curves correspond to the evolution for three labeled values of the pulse duration when 
F, = 20. 



In order to show the effect of ionizing pulses. Fig. (6) shows the evolution at constant pressure of a warm gas 
initially in thermo-chemical equilibrium close to the critical state for spontaneous condensation. The pressure Pq is 
bellow Pmax by a factor rj = Po/Pmax{Coo) = 0.8, when the pre-flash primary ionization rate is Coo = 10~^^s~^. The 
three curves in Fig. (6) correspond to the evolution for three different values of the pulse duration At = 0.7, 0.8, and 
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Fig. 7. The critical values of Fz as function of At for the three labeled values of r\. 

2.0 Myr, when the flash is initiated at t = 10^ yr and the factor of increment of C,o is = 20. Note in Fig. (6) that 
there is a critical value of At bellow which condensation does not occur (in this case ~ 0.7 Myr). Note also that the 
time required to complete the phase transition decreases as Aj increase, but the time required for the temperature 
drop from ~ 5000K to ~ 50K is insensitive to Aj. The kink at the bottom of the A* = 0.8 and 2.0 curves, and the 
smooth dip in the At = 0.7 curve are due to the inertia of the ionization fraction (the recombination time is much 
longer for the conditions corresponding to the upper curve). 

In Fig. (7) the critical values of Fz are plotted as function of At for the three labeled values of rj. These curves 
divide the plane (F^, At) in two regions: above the curves the flash is capable of inducing phase transition, and bellow 
the curve the gas returns to its initial state after the passing of the flash. Note that a cosmic ray pulse with a factor 
of increment F^ of the order of 10, and a duration of about one Myr can induce the condensation of the warm gas in 
TCE at 77 > 0.8. 

5. Summary and Conclusions 

The variation of the ionizing flux as a mechanism for stimulating the condensation of the diffuse gas was considered. To 

illustrate this effect, two situations were examined: one on the context of pregalactic conditions (a free metal cloud), 
and the other on the context of the actual interstellar medium (a gas with solar abundances). We have focused our 
attention on flash-like variations; that is, during a "short" period of time the ionizing flux is enhanced in comparison 
to the pre and post flash values. In both cases the cause of the induced phase change is the same: the enhancement of 
the cooling rate by the increase of the electron density caused by the momentary increase of ionizing flux. After the 
passing of the flash, the cooling rate remains enhanced due to the inertia of the ionization. In the first case (metal free 
gas) we show that after the lapse of UV flux decrement, recombination continues, but at a lower rate than the cooling. 
The excess of electrons at these relatively low temperatures results in an enhancement of the H2 rate formation due 
to the enhanced abundance of the H~ intermediary. Even when €2 > eo, the abundance of H2 may reach a large 
enough value to produce considerable self-shielding. If the cloud reaches a critical value (r ~ 1) for the optical depth 
at dissociating frequencies, the H2 abundance grows very fast, allowing the cooling of the cloud to temperatures of 
the order of lOOK. The temperature drop occurs in a fraction (~ 0.2) of the free fall time provoking a rapid decrease 
of the Jeans mass. However, if the post-flash dissociating level is large enough the H2 formation can be inhibit and 
the cloud remains warm. 

In the second case (solar abundances gas) the dominant cooling mechanism of the warm neutral gas (the excitation 
of heavy ions by electron impacts) is proportional to the electron density, and therefore, the ionizing flash increases 
the cooling efficiency. We considers fiash-like variations of the primary ionization rate by cosmic rays and calculate the 
marginal flash characteristics to induce warm gas condensation. We show that, for the expected states of the warm 
interstellar gas, ionizing flashes may easily induce the phase transition from the warm to the cool phase. The phase 
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transition is completed in about 10^ yr; however, the drop in the temperature from ~ 5000K to ^ 50K occurs in about 
5 X 10^ yr. 

The resuhs indicate that the mechanism of induced condensation studied here might play a relevant role in the 
gas evolution of the diffuse gas in both, the pregalactic and the actual interstellar medium conditions. 

The above results include only local processes. However, non-local processes like thermal conduction impose re- 
strictions on the size of the condensing structures. Thermal conduction tends to attenuate temperature gradients, and 
therefore, it imposes a minimal value for the mass of the condensing structures (Corbelli & Ferrara 1995; Ibahez & 
Rosenzweig 1995; Steele & Ibahez 1997). Critical masses of the order of galaxy masses are obtained in the case of a 
metal-free gas (Ibaiiez & Parravano 1983) and in the case of an actual interstellar gas, masses of the order of one solar 
mass are obtained (Parravano 1986; Parravano, Ibahez & Mendoza 1993). However, a wide range of critical masses 
are obtained depending on the initial gas state, and ambient conditions. Therefore, in addition to the restrictions on 
Fz and At to induce the warm gas condensation, the cosmic ray flash should cover a region greater than the critical 
value imposed by thermal conduction and other diffusive processes. 
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